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Abstract 

Using the general, model independent form of effective Hamiltonian, the decay 
spectrum of B° — > t:^ it~ l~~ decay is studied. The sensitivity of experimentally 
measurable asymmetries to the new Wilson coefficients and time is studied. It is 
observed that different asymmetries are sensitive to the new Wilson coefficients and 
they can serve as an efficient tool for establishing new physics beyond SM. 
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1 Introduction 



CP violation has been observed in the neutral kaon systems ]T|, Great effort is devoted 
to the study of possible signals of CP violation in B system, which will provide invaluable 
information about the origin of CP violation and is one of the most promising research area 
having the potential of establishing new physics beyond the standard model (SM). Started 
operating, two B-factories BaBar and Belle open an excited era for a comprehensive study 
of B meson physics and especially its rare decays. The main physics program of these 
factories constitutes a detailed study of CP violation in Bd meson and precise measurement 
of rare flavor-changing neutral current (FCNC) processes. Furthermore, the ultimate goal 
of these studies is to look for the inconsistencies within the SM (see for example H), in 
particular, to find indications for indications for new physics in the flavor and CP violating 
sectors. The factories mentioned above have both already signaled the first evidences about 
the CP violation in neutral B meson decays [f§|. 

One efficient way for detecting CP violation in the neutral B meson decays is the well 
known Dalitz plot asymmetry ||. The main advantage of the Dalitz plot asymmetries 
from the partial rate asymmetries is that they might be present even when partial rate 
asymmetries vanish. It should be noted that such CP-violating asymmetries in the angular 
variable ip (p is the angle between the £ + £~ and 7r + 7t~ planes) in the K — > tt + tt~£ + £~ 
decay was measured 0. In a recent report, NA48 Collaboration || confirms large CP- 
violating effect in the K — > tt + 7i^£ + £~ decay. Time evolution of the decay spectrum in 
K° (K° — ► 7r + 7r~ £ + £r) was investigated in framework of the SM in 0. 

In this connection several interesting questions come into mind immediately: how does 
the angular variable asymmetry A v evolve with time for the neutral B meson decay? How 
do new physics effects change the Dalitz plot asymmetry? Can we measure this asymmetry 
in 5-factories? This paper is devoted to answering these questions by analyzing the time 
dependence of the decays B°, B° — > ty + tt~£ + £~ in detail. It should be noted that the 
B — > n + n^£ + £^ decay in SM and beyond were studied in || and ||, respectively. 

The paper is organized as follows. In section 2 we present the general expression of 
the Dalitz plot asymmetries for the B° (B° — > 7r + 7T~£ + £~) decays using the most general, 
model independent form of the effective Hamiltonian. In section 3 we study the sensitivity 
of the Dalitz plot asymmetries on the new Wilson coefficients. 



2 Time evolution of the decay spectrum of B° — > tt + 7t £ + t 

The matrix element for the B — ► p{— > 7r + 7r~)£ + £~ decay is described by the b —>■ d£ + £~ 
transition at quark level. Following |j9|-[13|, the effective Hamiltonian of the b — > d£ + £~ 
transition can be written as the sum of the SM and new physics contribution, in a general 
model independent way in the following form 



Got \ - q v q v 

~l=-V td Va\ C S l dia^— Lb£^£ + C B r dia^— r, u > „ « 

V27T [ q q 
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+C LRRL d L b R £ R £ L + Crlrl d R b L £ R £ L + C T dcr^b £a^£ 
+tC TE e^ a(3 da IMU bia a ^ , 

where the chiral projection operators L and R in (|l|) are defined as 
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and Cx are the coefficients of the four-Fermi interactions. Here q = p + + p- and p + and p_ 
are the four momenta of £ + and £~, respectively. The first two of these coefficients, Csl and 
C BR: are the nonlocal Fermi interactions which correspond to -2m s C e 7 ff and -2m b C e 7 ff 
in the SM, respectively. The next four terms with coefficients Cn, Cl R , C R l and C RR 
describe the vector type interactions. Two of these vector interactions containing C L °[ and 
C L ° R do already exist in the SM in combinations of the form (Cg^ — Cio) and (Cg + Cio). 
Therefore, C L °l and C L ° R describe the sum of the contributions from SM and the new physics 
and they are defined as 



LL 
U LR 



C| + Cio + Clr ■ 



(2) 



The terms with coefficients Clrlr, C R llr, Clrrl and C R l R l describe the scalar type 
interactions. The remaining last two terms leaded by the coefficients CV and Cte, obviously, 
describe the tensor type interactions. 

The matrix element for the the process B — > 7t + 7t~£ + £~ can be obtained from the matrix 
element B — > p£ + £~, by the subsequent decay of the p meson to pion pair whose interaction 
Hamiltonian is 



(3) 



where is the the polarization four vector of the p meson, Q = p w + —p w - and p n + and p w - 
are the four momenta of 7r + and 7r~, respectively. Hence, in calculating the 4-body decay 
amplitude for the B — > tt + tt~£ + £~ process, the matrix element of the B — > p£ + £~ decay is 
needed. For this purpose it is enough to sandwich the effective Hamiltonian (|IJ) between 
the initial and final meson states. In other words, the following matrix elements 



^7m(! ± 75)& 



B 



d{l± l5 )b 



B 



derail ± 7 5 )& 



B 



need to be calculated in order to calculate the decay amplitude for the B — > pt 1 
These matrix elements can be written in terms of the form factors as 

p(P,e) d 7M (l± 75 )6 B(p B 



process. 



tti b + m p 



±ie*(m B + m p )Ai(g 2 



(4) 



Ti{PB + P\{e* q )^p- =F iq^{e*q) U 3 (g 2 ) 
m B + Tn p q A L 



Ao(q 2 ) 
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Multiplying both sides of Eq. (|J) with g M and using equation of motion, we get 
[p{P,e) d(l± l5 )b B(p B )\ 



(5) 



% 



[e*q){ =F (m fl + m p )A x (q 2 ) ± (m| - P 2 ^ M ' r) 



m B + m p 



±2m r 



A 3 (q 2 )-A (q 2 



The remaining matrix element is defined as 



da pv b 



B{ P i 



ie^J - 2T 1 (q 2 )e* x (p B + P) CT + -(m 2 B - P 2 ) T x (q 2 ) - T 2 (q 



e* X q a 



(6) 



T x {q 2 )-T 2 {q 2 )- 



m 2 B - P 2 



Uq 2 



(e*q)P x q a 



Moreover, the matrix element i^p(P,e) da^j^b B{p B )j can be calculated from Eq. 
using the identity 



°V75 — —-^iivafi® 



ai3 



which yields 

[p(P,e) da^jsb B(p B ] 

2{T 1 (q 2 )\e*(p B + P) u -e:(p B + P) IM - ^-{m% - P 2 )Kq u - e*^) T x {q 2 )-T 2 



+ - 



m 2 B — P 2 



Ts(q 2 



T x (q 2 )-T 2 (q 2 
Using Eqs. for the matrix element B — > -k + ti 



e*q)(P' M <f-P v q' M )j . 

£~ decay we get 
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where Vl s and V Ri are the coefficients of left- and right-handed leptonic currents with 
vector structure, respectively, Sl,r are the coefficients of the leptonic scalar currents and 
the last two terms correspond to the contribution of the leptonic tensor currents. These 
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new coefficients are functions of the Wilson coefficients and the form factors introduced in 
defining the hadronic matrix elements above. Their explicit expressions are given by 



Vl 2 
Vl 3 
Vl 4 
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I A, - C RL) — 



2mp (A 3 -A )-2(C BR -C SL )^ 





= v Ll {cr L - 


(~ltot 

U L R 1 


C R L 


— ► Crr) , 




Vr 2 


= v L2 {c L ot L - 


U L R J 


C R L 


- ► C RR ) , 




v Ri 


= vuczi- 


/~itot 
~* °L_R J 


C R L 


c RR ) , 
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where dependence on q 2 is implied. To obtain the amplitudes for the B° 
decay, it is enough to use CPT invariance and make the replacements in Eq 



2,3,4) 



- 2m p (A 3 - Aq) 

- 2m p (A 3 - Aq) 

1T + 1T~£ + 1 





-> —Vlv,Ri j 


Vhi;Ri 


■* V Li;Ri , (i 


SL; R ~ 


Sr,L , 


Bqj-s - 


^6;7;8 , 


-Bg;10;ll ~ 


~Bg-iQ. n , 
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where the meaning of the bar can be explained as follows: If any arbitrary one of the 
coefficients Vl^iu * s represented in the form 

Vl iA = \V Li ., Ri \e i5 ^e i<f,L ^ , 

where 5 and are the strong and the weak phases, respectively, then 

V Ll - Rl = \V Li , Ri \e lS ^e- 1 ^ . 

The starting point in analysis of the time evolution spectrum of the B°(B°) 
decay is Eq. (H). The time evolution of the B mesons is given by 



7T 1 7T 



B°(t) 
B°(t) 



9+(t) 
9+(t) 



B° 
B° 



&-(*) 



B° 
B° 



with 



g± = exp 



B 



— irriB]t 



cos 



Am 



/Am 



sm 



I" 



-t 



and £ = p/q. Here m^(r^) and Am^Ar^) are the average and the difference of the 
masses (widths) of the two mass eigenstates Bh and B^, respectively. 



The matrix element of the time dependent B°(t) 
from Eq. (Rf) by making the following replacements. 



7T 7T 



decay can be obtained 



^9:10:11 ~~ > 



>v Ll ., Rl (t) 

*V Li]Ri {t) 
B 6 . 7; s(t) 

Bg-I0;ll(t) 



Vl,;^ + -zg-(t)V Li;Ri , 
g+(t) S L . R + -g„(t)~5 R . !L , 

g+(t) 5 6;7;8 + | i?6;7;8 , 
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9+{t)B, 



9:10:11 



7 9-(t) ^9;10;ll 



Similarly the matrix element for the B°(t) 
respective matrix element Ai(B° 



ir+ir-e+t 



7T + 7r £ + £ decay can be derived from the 
) with the help of the replacements 



— > BQ- 7 -s{t) 



B. 



9;10;11 



^9;10;11 if) 



g+{t) V Ll - Rl -£g-(t) V Ll - Rl , 
2+(t)V L!;R! +^_(t)V LlA , 
g+(t)S L . R + £g_(t)S R . L , 
g+(t) B 6 . 7 . s + £ g_(t) B 6 . 7 . s , 
g+(t) B g . 10 . n - £ g_(t) B 9 . 10 . n 



2,3,4) 
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Using the formalism for the Ki A decay ||14[| , from the matrix element (§) one can obtain 
the B° — > 7t + tt~£ + £~ decay rate in terms of the following five variables: invariant mass 
Sm = P 2 = (p n + + Ptt-) 2 °f 7T + 7r _ pair, invariant mass Se = q 2 = (p + + pJ) 2 of £ + £~ pair; 
the angle 9 between p n + and p + + measured with respect to the center of mass of 7r + 7r~ 
pair; the angle 9g between p + and p w + +p n - , measured with respect to the center of mass of 
£ + £~ pair and ip is the angle between the normals to the 7r + 7r~ and £ + £~ planes. The final 
four-body phase volume in terms of the above-mentioned five variables can be written as 



dX 



PS 



2(47r) 6 m| 



XP/3e dsM dsg d(cos 9) d[cos 9g) dip 



(12) 



where 



X = 


(Pq) 2 - s M s L 


(3 = 


— \ 1/2 (s M ,ml 
sm 


fit = 


— X 1/2 (si,mj,m 


St 



1/2 
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The bounds of the integration variables are 



4774 < s m < m 2 B , 

Am] < s e < (m B - y/slf) 2 , 

o < e, Bi < 7T , 

< ip < 2tt . 

Different scalar quantities appearing in expression of the |.M| 2 can be expressed in terms 
of the above-mentioned five variables as 



Pq 


= \( m l -s M - s e ) , 


QN 


= (3f3e (Pq cos 9 cosBg — ^s M si sin 6* sm9e cosp>) 


PbP± 


= ^(s M + Pq±X(3 cos6) , 


qN 


= o, 


Qq 


= XpcosO, 


QP 


= o, 


NP 


= X [3 e cos e , 


Q 2 


= smP 2 , 


N 2 


= sift* , 


e»vX„Q»P u NY 


= —^smSi X [3 Pi sin 9 sin 9e sin p> , 



where 



N = p + — p_ 



The differential decay rate for the B°(t) — > tt + tt~£ + £~ is (in the expression below leptons are 
taken to be massless. However in calculation of the decay spectrum for B°(t) — > 7r + 7r~T + T~, 
the r lepton mass is taken into account. But this expression is quite lengthy and for this 
reason we do not present it in the text) 



dY = 



Ga 



9, 



piT1T 



(sM-mlf + rlml 2 14 7T 6 777 3 B 



4V2n tbVtd 
xdsM dsed(cos9) d(cos9e) dip X , 

where 

X — Xi cos 2 9 cos 2 9e + X 2 cos 2 9 cos 9i + J 3 cos 2 # + T 
+J 5 sin(2#) sin(20*) cos^ + J 6 sin(2#) sin(2^ 
+X 8 sin(26>) sin^ cosy? + X 9 cos 2 9 e + Z 10 cost 



X/3 



(13) 



and 



sin 2 6> sin 2 ^ sin(2</?) 
h) sin(/9 + X 7 sin(26>) sin 6^ simp (14) 
(p + X 9 cos 2 9 e + J 10 cos 0^ + Jn sin 2 # sin 2 ^ cos 2 y9 + J 1; 
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x 2 
1* 



X x = (3 2 {x 2 [\V Ll \ 2 -\V L3 \ 2 + \V Rl \ 2 -\V Rs \ 2 ] -2X(Pq)[2Re(v L2 Vl s 
+ 2Re(v R2 V Rs )} -4(Pg) 2 [A(|V il | 2 + |V i?1 | 2 ) + |V i2 | 2 + |V i?2 | 2 ]} , 

2f3 2 V\[\ - 4 [Pqf] [ - Re (V^V^) + Re (V Rl V R2 ) ' , 

(3 2 \{\ (|V L3 | 2 + \V R f) + |V fi2 | 2 - 4 (Pg) [Re (v L2 V2 3 ) + Re (Vr 2 V R3 ) ] } , 

4/3 2 v / A[lm (V L2 V£ 1 ) + Im (Vr 2 V* Li ] 



A 



Re (V L2 V1 3 ) + Re (V R2 V Rs ) ] - 2 (Pq) (|V La | 2 + \V R f) } 



= (3 2 V\{\[lm (V L3 V* Ll ) + Im (V^V^) ] + 2 (Pg) [im (V Zl V£ 2 
+ Im(v fll V* 2 )]}, 

= 2(3 2 X v^i^ [im (V L3 V1 2 ) ~ Im (Vr,V* R2 



= 2(3 2 \J\~sms~£ I A [Re (v Ll V2 3 ) - Re (v fll V^ 3 ) ] + 2 (Pg) [ - Re (V Ll V* L2 
+ Re(v fil V^ 2 )]}, 



J 9 = 4/? 2 ^ M A[|V Ll | 2 + |V /?1 | 2 



4p 2 s e s M A (|V Ll | 2 + |V Rl | 2 ) - \V L2 \ 2 - \V 



y R 2 \ 



Z10 
In 

X 12 = 4(3 2 s e s M [\V L2 \ 2 + \V R2 
Taking the narrow resonance limit of p meson, i.e., using the equations 



lim 

r- 



9%*™ p f Am 2 \ 3/2 
48tt V m 2 ' 



and, 



(s M - my) 



+ ° (s M - my) + m 2 p 

we perform integration over sm easily in Eq. ([Hp. After integrating over sm, si and 9, we 
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get for the angular distribution of the B°(t) — > 7r + 7r £ + £ decay 

2 



d(cos8i) dip 



Ga 



4v/2' 



7T 



v tb v; d 



48vr 2 



4m^ / m 2 , 



1/2 



x 



| (Jx{ + 2X^ cos 2 ^ + (Jx^ + 2X^) cos0< + ~X 4 sin 2 ^ sin(2<^) 



+ sin 2 ^ [l + cos(2<^) 
3 



; x' + 2x: 
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(15) 



where we have introduced the notation X- = J Xlidsg. CP violation manifests itself in the 
coefficients ~ 2 4 and ~ X' u . The dominant CP violation term is contained in the coefficient 
(4/3)24. After integrating over cos 6^ we obtain the 92-distribution 



dT 

dip 



2 /2 



3 V3 



16 



Z[ + 21' 9 ) + — X 4 sin(2p) + -Z' n 1 + cos(2<^) + 2 [-1' 3 + 2X( 2 



9 
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As has already been noted the term ~ sin(2<^) is CP- and T-odd. Hence this term can be 
isolated by considering the following asymmetry in terms of the differential decay rates of 
the B meson with respect to ip 



A v (t) 



The main feature of such CP violating asymmetries is that they can be obtained by con- 
sidering the sum of the differential decay rates of the B and B rather than the difference 
of these rates, as usual. 

Furthermore we study the dependence of the asymmetry between the spectrum inte- 
grated rates of B°(t) — > tt + tt~£ + £~ and B°(t) — > tt + tt~£ + £~ decays, on the new Wilson 
coefficients, which is defined as 





/■tt/2 

1 -J 


r +, 

tt/2 J 


rZir/2 

l 'J 


f" 1 

3tt/2_ 




I, 


«r/2 


r +, 

tt/2 J 


p3tt/2 

, + J 


f 2, ' 
3n/2^ 





A^t) 



T(B°(t) 



7T 7T 



T{B°(t) -> 7T + 7T- 



(16) 



T(B°(t) -> 7T + 7T-£ + r) + r(B°(t) -> 7r + 7rl + ' 

It should noted that it is also possible to construct a differential asymmetry that isolates 



another combination of the terms containing imaginary parts (see for example [15, [If 
Such terms, for example, can be isolated by considering the difference distribution of the 
same hemisphere and opposite hemisphere events. The corresponding asymmetry can be 
defined then as follows 



A 2 {t) 



1 

r 



2- 



dip 





-1 



d(cos9g 



+1 



G?(cOS 9) 



dT 



X 



dip d(cos 61) d(cos 9) 



(17) 
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3 Numerical analysis 



In this section we study the dependence of the asymmetries A v , Ai and A2 on new Wilson 
coefficients and time for the B°(t) — > 7t + tt~£ + £~ decay. Before presenting the numerical 
results we want to note that, although the expressions for the asymmetries A v , A\ and A2 
are given for the most general case, i.e., each one of new Wilson coefficients might have new 
strong and weak phases. In our numerical analysis we choose all new Wilson coefficients to 
be real and they all vary in the region between —4 and +4. 

For the form factors which are needed in the course of performing numerical calculations, 
we have used the prediction of the light cone QCD sum rules. In our numerical analysis we 
will use the results of the work |17j (see also [[18], [19]]) in which the form factors are described 
by a three-parameter fit where the radiative corrections up to leading twist contribution 
and SU(3)-breaking effects are taken into account. The g 2 -dependence of the form factors 
which appear in our analysis could be parametrized as 



F(s) 



F(0) 



1 — dp s + bp s 2 



where s = q 2 /m 2 B is the dilepton invariant mass in units of B meson mass, and the param- 
eters F(0), dp and hp are listed in Table 1 for each form factor. 





F(0) 


dp 


bp 






0.372 


1.40 


0.437 


aT 


^p 


0.261 


0.29 


-0.415 


At 




0.223 


0.93 


-0.092 


yB- 




0.338 


1.37 


0.315 


T B- 


->p 


0.143 


1.41 


0.361 


Tf~ 


^p 


0.143 


0.28 


-0.500 




*p 


0.101 


1.06 


-0.076 



Table 1: The form factors for B — > p£ + £ in a three-parameter fit. 
As it comes to the values of the Wilson coefficients C7 \rrib) and Cio(rrib), whose an- 



alytical expressions in the standard model are given in [5D, 2I|, are strictly real as can be 
read off from Table 2. In the leading logarithmic approximation, at the scale 0(fi = mj), 
we have 



C e 7 ff (m b ) 



-0.313 
-4.669 



Although individual Wilson coefficients at fi ~ m,b level are all real (see Table 2), the 
effective Wilson coefficient Cg^ (nib, s) has a finite phase, and in next-to-leading order 
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Ci 


c 2 


c 3 


c 4 


c 5 


c 6 


c 7 e// 


c 9 


r< e ff 


C(o) 


-0.248 


1.107 


0.011 


-0.026 


0.007 


-0.031 


-0.313 


4.344 


-4.669 


0.362 



Table 2: The numerical values of the Wilson coefficients at \i ~ m;, scale within the SM. 



C e 9 ff i 



m b , s) 



C 9 (m b ) 



a 8 (fi) 
1 H w s 



7T 



+ s) + Y LD (m B , s) 



where C 9 (m b ) = 4.344. Here u (s) represents the 0(a s ) corrections coming from one-gluon 
exchange in the matrix element of the corresponding operator, whose explicit form can be 



20, El]. Here Y$d can be obtained by a 



found in P0fl . In (|TJ) Y SD and Y LD represent, respectively, the short- and long-distance 
contributions of the four-quark operators O i= i t . 
perturbative calculation 

Y SD (m b ,s) = g(m c ,s)C^-^g(l,s)[4C 3 + 4C A + 3C 5 + C 6 ] 
- lg(p,§)[C 3 + 3C4\ + ~[3C 3 + C A + 3C s + C 6 ] 



- X u [3d + C 2 ] [g (0, s) - g (m c , s)) , 



where 



C (0) 



3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 



v tb v; d 



and the loop function g (m q , s) stands for the loops of quarks with mass m q at the dilepton 
invariant mass s. This function develops absorptive parts for dilepton energies s = Am 2 : 



9 {m q , s) = -- lnm g + — + -y q - - (2 + y q ) J|l 



1) 2 arctan 



y/Vq ~ r . 



where m g = m q /mb and y 9 = 4m^/s. In addition to these perturbative contributions 
cc loops can excite low-lying charmonium states ip(ls), ■ • • , tp(6s) whose contributions are 
represented by Yld [P2"| . However in the present work we restrict ourselves to the consid- 
eration of short distance contributions only and therefore we find it redundant to display 
the explicit form of Yld- 

Calculations show that the asymmetry A v for the B°(t) — > 7r + 7r~~/i + //~ decay is sensitive 
only to tensor type interaction. We observe that, practically, the existence of the new 
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vector or scalar type interactions do not effect the SM results crucially, which is about 
13%. However when tensor type interactions are taken into consideration this asymmetry 
is reduced noticeably and (Ap) ~ 2.5%. Therefore if in future experiments the results 
obtained results differ from the SM results, it is an unambiguous indication of the existence 
of new physics beyond SM. 

As one considers the B°(t) — ► 7t + 7t _ t + t~ decay, the maximum value of the asymmetry 
A v is about 5%. Similarly, as in the B°(t) — ► it + n~ n + [i~ decay, A^ is not sensitive to the 
new vector and scalar interactions and for the contribution of the tensor interaction is less 
1%. 

We further studied the dependence of the asymmetry A\ on the new Wilson coefficients 
and time for both B°(t) — > tt + tt~ fi + fi~~ and B°(t) — > 7t + tt~t + t~ decays. Our calculations 
show that the maximal possible value that A% asymmetry gets for the B°(t) — > 7r + 7r~/i + /i~ 
decay is about 2% for all type of new interactions. Similar behavior takes place for the 
B°(t) — > 7t + tt~t + t~ case as well, with only one exception. For tensor type interactions 
which are controlled by the new Wilson coefficients CV and Cte, A\ can arrive at large 
values about 20% (see Figs. (1) and (2)). So any departure from the SM prediction in 
future experiments, is an indication of the presence of tensor type interactions. 

Finally we have studied the dependence of the Ai asymmetry on new Wilson coefficients 
and time for the above mentioned decays. For the B°(t) — > tt + tt~ fi + fi~ decay the maximal 
value of Ai is about coefficient Cll = —4 causes Ai to decrease about 50% from the SM 
prediction, but if Cll = +4, practically the value of A2 seems to coincide with the SM 
prediction, which also gives us a clue in determining the sign of of the new vector interaction, 
as can be seen in Fig. (3). As is depicted in Fig. (4), similar situation seems to hold for 
the new Wilson coefficient Cm- Our numerical analysis shows that A2 is quite insensitive 
to the existence of scalar interaction and yields almost same results with SM prediction. 
But an important observation is that the tensor interaction enhances the SM prediction by 
about 50%. Therefore it would be an unambiguous confirmation of the existence of new 
vector (tensor) interaction, if in future experiments the measured value of Ai appears to 
be less (larger) than the SM prediction. 

On the other hand for all values of the new Wilson coefficients, the A% asymmetry for 
the B°(t) — > 7t + 7t~t + t~ case coincides almost with the SM prediction. For this reason A2 
for B (t) — > 7t + tt~t + t~ decay does not seem to provide us with any new information about 
new physics beyond SM. 

As a final remark we should emphasize that similar analysis can be performed for the 
B°(t) — > Tr + 7Y~£ + £~ decay as well. For this aim it is enough to use Eqs. (§), @ and ([TT]) • 



In conclusion, in this work we have studied the time evolution of the decay spectrum 
for the B°(t) — ► tt + tt~£ + £~ decay. The sensitivity of the experimentally measurable asym- 
metries Ap, Ai and A2 to the new Wilson coefficients and time is studied in detail. It is 
observed that different asymmetries show different dependencies on different new Wilson 
coefficients for the B°(t) — > 7t + tt~ fi + and B°(t) — ► 7r + 7r~T + T~ channels. Studying differ- 
ent asymmetries for the above-mentioned two channels on new Wilson coefficients can give 
essential information about new physics and can serve as an efficient tool in determining 
not only of the magnitude, but also of the sign of new Wilson coefficients. 
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Figure captions 

Fig. (1) The dependence of the asymmetry A\ on the new Wilson coefficient Ct and time 
for B°(t) — > it + it~t + t~ decay, where t s is the B° meson life time. 

Fig. (2) The same as in Fig. (1), but for the new Wilson coefficient Cte- 

Fig. (3) The dependence of the asymmetry A2 on the new Wilson coefficient Cll and 
time for B°(t) — > ir + iT~fi + fi~ decay. 

Fig. (4) The same as in Fig. (3), but for the new Wilson coefficient C RL . 
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